#ifndef LINEARSYSTEM_H
#define LINEARSYSTEM_H

#include <boost/numeric/ublas/io.hpp>
#include <boost/numeric/ublas/lu.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/vector.hpp>

typedef boost::numeric::ublas::compressed_matrix<double> compressed_matrix;
typedef boost::numeric::ublas::compressed_vector<double> compressed_vector;

void solveLinearSystem(compressed_vector& x, compressed_matrix& A, compressed_vector& b);

#endif
